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ABSTRACT 

We present measurements of relative elemental abundances in plumes and inter¬ 
plumes. Plumes are bright, narrow structures in coronal holes that extend along open 
magnetic field lines far out into the corona. Previous work has found that in some coro¬ 
nal structures the abundances of elements with a low first ionization potential (PIP) 

< 10 eV are enhanced relative to their photospheric abundances. This coronal-to- 
photospheric abundance ratio, commonly called the FIP bias, is typically 1 for element 
with a high-FIP (> 10 eV). We have used EIS spectroscopic observations made on 2007 
March 13 and 14 over an ~ 24 hour period to characterize abundance variations in 
plumes and interplumes. To assess their elemental composition, we have used a differ¬ 
ential emission measure (DEM) analysis, which accounts for the thermal structure of 
the observed plasma. We have used lines from ions of iron, silicon, and sulfur. Erom 
these we have estimated the ratio of the iron and silicon EIP bias relative to that for sul¬ 
fur. From the results, we have created FlP-bias-ratio maps. We find that the FlP-bias 
ratio is sometimes higher in plumes than in interplumes and that this enhancement can 
be time dependent. These results may help to identify whether plumes or interplumes 
contribute to the fast solar wind observed in situ and may also provides constraints on 
the formation and heating mechanisms of plumes. 

Subject headings: Sun: corona - Sun: UV radiation - Sun: elemental abundances - solar 
wind - plumes 


1. Introduction 


Polar plumes are bright, narrow structures rooted in coronal holes. They extend along open 
magnetic field lines, expanding with alti tude as the magnetic held dive rges. In white light, plumes 
have been observed as radial structures (|van de HulstlllOSQ : ISaitolll965ll that reach altitudes greater 
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than ~ 3Oi?0 (jPeForest et alJl200lh . Plumes are generally denser and cooler than their surround¬ 
ings and are likel y heated at the base (IWandll994l i . An extensive review of plume plasma properties 
can be found in Wilhelm et al. I (HT Despite many studies carried out in white light and ex¬ 
treme ultraviolet (EUV), there are still many open questions about plumes, such as their formation 
mechanism and their possible contribution to the fast solar wind. 

Because of projection effects, the three-dimensional nature of plumes is still controversial. 
The classical view considers plumes to be quasi-cylindrical structures formed by a small emerg¬ 
ing magnetic bipole reconnecting with surrounding open magnetic field lines. Such plumes are 
kn own as “beam ” plum es. A second type of plume, called a “curtain-plume”, has been proposed 
by iGabriel et al.l (|2003l i . These plumes correspond to a chance alignment along the line of sight 
(LOS) of a series of the roughly vertical, faint outflows along the boundaries of the magnetic net¬ 
work on the solar surface. Since these two plume populations present almost the same plasma 


properties (|Gabriel et al.ll2009lL it is difficult to observationally distinguish the two types. 


It is now accepted that coronal holes are the origin of the fast solar wind ([Krieger et al.lll973l L 
However, it is unknown whether plumes or interplumes are the main contributor of material to the 
fast solar wind. Some recent results have found that the outflow velocity is smaller in plumes than 


in the surrounding interplume gas (jTeriaca et al. 


2003 


Raouafi et al.l 120071 ). Since the solar wind 


speed at low heights is higher in interplumes than in plumes, these authors suggest that interplume 
regions are the preferred channel for the acceleration of the material into the fast solar wind. 
Gonversely, other studies measured higher outflow velocities in plume s than in interplumes ari d 
concluded that plumes are the main contributor to the fast solar wind, ( Gabriel et al. 20^. 2QQ^. 
The oretical mode l s of p lumes extending into t he solar wind have bee n developed bvlWand (|l994l ) 
and lHabbal et al.l (jl995l j and more recently bv iPinto et al.l (|2009l i and IVelli et al.l (120111 1. However, 
neither theory nor observations have been able to conclusively determine the main contributor to 
the fast solar wind. 


One way to determine if plumes or interplumes feed the fast solar wind is by searching for 
corresponding material in the in-situ measurements. Goronal abundances have been found to vary 
signihcantly among some coronal structures. Since the composition of the solar corona becomes 
hxed at low heights, elemental abundances measurements at these heights can thus be potentially 
used to identify the sources of the fast solar wind. 


These coronal abundance variations are correlated with the first ionization potential (FIP) 
specihc to each element. The abundances of low-FIP elements (FIP < 10 eV) appear to be 
enhanced relative to their photospheric abundances in some coronal structures. On the other 
hand, elements having a high-FIP value (FIP > 10 eV) exhibit coronal abundances close to photo- 
spheri c. For an exhaus t ive re view of the abundance measureme nts in the solar corona see iRavmond 
(I 1999 I I , IFeldman et al.l (j2005l i , and iFeldman &: Widind (j2007l l . It is convenient to report this so- 
called FIP effect in terms of a FIP bias, which is dehned as 

/X — ; 
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where and are the coronal and photospheric abundances of element X, respectively. 


_ It is st ill unclear wh e ther o r not plumes are subject to a FIP effect. Measurements bv lWilhehn &: Bodmer 

( IQQsI l and Young et al. (1999) with the Coronal Diagnostic Spectrometer (CDS) on board the So¬ 
lar and Heliosph eric Observa t ory (S OHO) concluded that there is a FIP effect in plumes. A more 
recent study bv ICurdt et al.l (120081 ) using measurements from the SOHO/Solax Ultraviolet Mea¬ 
surements of Emitted Radiation (SUME R) spectrometer f o und a small enhancement in plumes. 
Conversely, using data from SOHO/CO'S. iDel Zanna et al.l (j2003j) found no significant EIP effect 
in plumes. A time dependence of the FIP bias in plumes could possibly explain the discrepancies 
among these studies. Such behavior h as been seen in measurem ents in active regions where the 
FIP bias gradually increased over time ((Feldman &: Widindl2003l b 


Here, we investigate the elemental composition of both plumes and interplumes and their 
evolution over an ~ 24 hour period. Our measurements are base d on a differential emission measure 
(DEM) analysis ( Eeldman &: Widingl (20031 : (Feldman et al.ll2009l b The DEM provides a measure of 
the amount of LOS emitting material as a function of electron temperature. Our DEM analysis 
was carried out using many emission lines from the low-FIP elements iron (EIP ~ 7.90 eV) and 
silicon (FIP ~ 8.15 eV). Relative abundances measurements are derived using the intermediate-FIP 
element sulfur (FIP 10.36 eV). 


The rest of this paper is organized as follows. In Section [2l we describe the observations and 
the data processing. The DEM analysis is reviewed in detail in Section [3l The results are described 
in Section [H Systematic uncertainties are reviewed in Section [5l The results are discussed in 
Section [6l Finally, we summarize in Section [71 


2. Data Analysis 

2.1. Plume and interplume observations 

We used off-lim b observatio ns from the Extre me Ultraviolet Imaging Spectrometer (EIS; 
Culhane et al.l l2007l ) on Hinode ((Kosugi et al.ll2007l ). The data were collected with the 1" slit. 
We used 6 observations covering the full wavelength range and spanning an ~ 24 hour period. The 
first observation was performed on 2007 March 13 starting at 12:00 UT and ending at 15:39. The slit 
was rastered from X = 130.57" to —124.13", respectively, and centered vertically at Y = —972.43", 
with an exposure time of 50 s for each of the 256 pointings. Eigure[T] displays the positions of the slit, 
overlaid on the closest 171 A image from the S'ORO/Extreme Ultraviolet Imaging Telescope (EIT; 
Delaboudiniere et al.lll995l L We also used four smaller rasters made on this same day, starting 
respectively at 19:05, 20:07, 21:08 and 22:10. For these data, the 1" slit was centered vertically at 
Y = —972.42" and scanned from X = 30.80" to —28.01" with a 60 s exposure time for each of the 
60 pointings. The last observation was performed on 2007 March 14 from 11:05 to 14:44, with the 
same properties as the first observation. These off-limb rasters are ideal for our study because they 
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present a well developed polar coronal hole and contain insignificant contamination from quiet Sun 
plasma (see Section 13.41 for more details). 


The data were processed using standard routines available in the Interactive Date Language 
(IDL) Solar Software (SSW) package. CCD dark current, spikes, and warm pixels were removed and 
the intensity was calibrated wit h the eis_prep . pro routine. Next, missing pixels were interpolated, 
following the implementation of Young ( 201ol i and using the eis_replace_missing.pro procedure. 
Systematic drifts in the waveleng th scale were t hen t aken into account using eis_shift_spec.pro, 
based on the empirical model of iKamio et al.l (|2010|). The four smaller rasters on 2007 March 13 
were coaligned and averaged together in order to improve their statistical accuracy. Finally, each of 
the three resulting data sets was summed in both the horizontal and vertical directions into spatial 
bins of 4 pixels by 4 pixels. 


The fields of view seen by FIT in the 171 A bandpass corresponding to our three data sets 
are shown in the left column of Figure [2j In the first observation, two different plumes are clearly 
visible in the center of the raster, another one on the right hand side, and a suggestion of a fourth 
one is visible on the left. One plume seems to come from a bright point, whereas the origin of the 
others is difficult to determine. The field of view is greatly reduced in the second data set, but one 
plume still falls within the raster. In the last observation, only three plumes are observed; the one 
coming from the bright point has completely disappeared. 


All of the observations used the full spectral range of EIS. This, combined with long expo¬ 
sure times, provided usable spectral lines covering a wide range of electron temperatures, between 
logTe = 5.7 to 6.3. Here and throughout, temperatures are given in units of Kelvin. We used 4 
lines from Si VII, 2 from Si X, and 15 from Fe charge states spanning VIII-XV. These lines are 
mostly unblended (see Tabled]). Measurements of relative abundances were derived using these 21 
lines combined with the S VIII 198.55 A line. 


2.2. Line fitting 


We ht Gaussian prohles to each spectral line in order t o derive inten sities. For this, we 
used the eis_auto_f it .pro procedure, described in detail bv lYound (j2013l l. Most of the lines 
chosen are unblended; in this case a single Gaussian fit was applied to the data. The initial 
EIS intensity error bars are derived from the IDL routine e is_prep.pro , whic h include photon 
counting statistics and error estimates for the dark current. iHahn et al.l (|2012l l have found that 
this approach underestimates the uncertainties and so we adopted their Monte-Garlo procedure to 
derive the uncertainties for each of the Gaussian parameters used in the analysis here. Note that 
here and throughout, all uncertainties are given at an estimated la statistical accuracy. 

The few blended lines in our data set were easily dealt with. The Si VII 278.45 A line is 
blended with the Mg VII 278.3 9 A line. The lines we re separated using a 2-Gaussian fit that forced 
them to have the same width ([Landi &: Youndl2009l L This assumption is reasonable because both 
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the corresponding thermal and non-thermal widths of each ion are expected to be similar. The 
Fe XII 195.12 A is self-blended with the Fe XII 192.18 A line, leading to a distortion of the p rofile 
from a single Gaussian. For this line, we follow the procedure described bv I Young et ahl (j2009|) and 
use a double Gaussian fit in order to derive the Fe XII 195.12 A intensity. The S VIII 198.55 A 
line is blended with the Fe XI 198.54 A line. However, the Fe XI contribution c an be estimated 
throu gh the theoretical Fe XI 198.54 A/189.12 A branching ratio, which is ~ 0.8 ( Landi fc Young 
20091 ). The Fe XI 189.12 A is unblended and relatively isolated from other lines, so the intensity 
can be accurately determined using a single Gaussian fit. We then estimate and subtract the 
Fe XI 198.54 A contribution to find the S VIII 198.55 A intensity. Fortunately, the Fe XI 198.54 A 
intensity is very weak, with a ratio of Fe XI/S VIII intensities which never exceeds 15% for off-limb 
observations. Taking into account the estimated 10% accuracy of the branching ratio, the resulting 
uncertainty in the final S VIII 198.55 A intensity does not exceed 2%. 

We also took into account instrumental scattered light in our analysis. This light is superim¬ 
posed on the true coronal line emission and may lead to an overestimation of the intensity, especially 
at high heights above the limb. To remove the stray light from our data, we subtracted a stray light 
line profile. The u pper limit of the off-limb scattered li ght is estimated to be 2% of the averaged 
on-disk intensities (jUgarte-Urral l201Cl : iHahn et al.ll2012l h Thus, we hrst fit the lines with a single 
Gaussian for on-disk data with R ^ 0.95 Rq, as measured from the Sun center. Then, we subtract 
from each off-limb spatial bin a Gaussian stray light line prohle using the same centroid position 
and width and 2% of the averaged on-disk intensity. The stray light in our analysed data does not 
exceed 2.5% of the total measured for any given line intensity. 


2.3. Density diagnostic 


Some of the lines used in this study are density sensitive, such as the Fe IX 189.94 A and 
Fe XII 195.12 A lines. The intensity for these transitions depend on both te mperature and density . 
This is a potentia l source of systematic uncertainty in the DEM inversion (jGraig &: Brownlll976l : 
Brown et al. 1991 : Judge et al. 1997 : Guennou et al. 2012al f5l. To mitigate this problem, we have 
deter mined the density using the intensity ratio of the Fe IX 188.50 A and 189.94 A lines (Young 
2009|). For this, the plasma is assumed to be isothermal, at a temperature about logTe « 5.9, 
corresponding to the peak of the contribution function of the two lines. We used the atomic data 
from GHIANTI version 7.1 ( Dere et'aLlll997 : Landi et al. 2013 ). The right column of Figure [2] shows 
the density measured for our data sets as a function of the solar-Y and Y coordinates. The black 
contours correspond to the intensity observed simultaneously in the Fe X 184.54 A line. The inferred 
density is 1.3 ± 0.2 x 10® cm“® in plumes and 7.9 ± 0.8 x 10^ cm~® in interplumes a t 77 = 1.02 77, 


a 


This agrees with typical values previously measured in plumes and interplumes (IWilhelm et al 

201lh . 
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3. Measuring FlP-bias ratios 


A DEM analysis, by accounting for the temperature dependence of the observed intensities, is _ 

one w ay to measure spectroscopically the elemental composition of an emitting plasma (jFeldman fc Widine 
2 OO 3 I I. For the work here, we first d erived the DEM usi n g the Fe lines. The DEM was reconstructed 
using the regularization method of lHannah fc Kontad (120121 ). Taking this DEM we predicted the 
Si line intensities and found good agreement with those observed. This implies that both ions 
have a similar FIP bias. Hence, we combined the Fe and Si data to compute the DEM for each 
spatial bin of our three data sets. To maximize the reliability of the DEM inversion, we imposed 
various criteria pertaining to the number, intensity, temperature coverage, and uncertainty of the 
spectral lines used (see Section EH). Based on this Fe & Si DEM, we then computed the predicted 
intensities of the moderate-FIP S VIII 198.55 A line. Finally, by comparing both the predicted and 
observed S VIII intensities, we determined the ratio for the FIP bias of the low-FIP elements Fe 
and Si relative to that for S. 


3.1. Differential Emission Measure (DEM) 

In an optically thin plasma, the intensity of a given spectral line corresponding to the transition 
from level j to level i for element X is given by 

Iji = ^ I G,i{Te,n,)nldl, ( 2 ) 

where dl is the differential along the LOS; Tg and Ug are the electron temperature and density, 
respectively; Ax is the abundance of X relative to hydrogen; and Gji is the contribution function, 
which is defined as 


Gji^ei'^e) — Aj 


nAX+™)n(X+" 


n{H) 1 


n(X+™') n(X) Ug Ug 


(3) 


Here, Aji is the radiative transition rate, nj(X'*“'”) represents the number density of ion X'*'™’ i 


m 


level j, n(X“*'™’) is the number density of ion X"*""*, n(X) is the number density of the element X, 
and n{H)/ne is the abundance of hydrogen relative to that of electrons. 


Following iPottaschl (jl963l . Il964l l , Equation ([2]) can be reformulated as 

r *+00 


Iji — 


^ r 

4vr Jo 


Gji{Te,ne)({Te)dlogTe, 


(4) 


where the DEM ^ provides a measure of the amount of emitt ing plasma as a functi on of the electron 
temperature. The DEM is defined as ^{Te) = nldl/dlogTe I Craig fc Brown 1976). 

Reliably inferring the DEM from a series of observations has proved to be a challenge, due 
to the inverse nature of the problem (|Craig &: BrownI Il976l : iBrown et al.l Il99ll : iTesta et al.l l2012l : 
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Guennou et al.ll2013l ). To ensure an accurate solution, many lines having low uncertainties and 


spanning a broad temperature range are needed. In order to automate the DEM inversion process 
for the large number of spatial bins in our three data sets, we have imposed several criteria to 
reject spatial bins for which the DEM is poorly constrained. We required at least 7 different 
lines to perform the DEM inversion for each spatial bin. In order to ensure adequate temperature 
coverage, we also require that there be at least 5 different charge states in the available lines used 
to compute the DEM. Since the EIP bias is deduced using the S VIII 198.55 A line, which forms 
at about log Eg = 5.8, we required at least 3 different charge state spanning the temperature range 
log Eg = [5.7,6.0]. Additionally, for each spatial bin, lines are rejected if they are too weak, which 
we define as Iji < 10 ergs cm“^ s“^ sr“^. Lines were also considered too noisy and rejected if the 
signal-to-noise ratio was less than 3, i.e., Iji < 3cj, where a is the uncertainty of Iji (see Section r2.2l) . 
All together these conditions ensure that the analysis uses only spatial bins for which the DEM is 
well constrained. 


3.2. The iron DEM 


The first step in our approach is to derive the DEM in the coronal hole using only the iron 
lines. The coronal DEM is related to the observed Ee intensities through 

AC r+oo 

i!L = 1? n,)UTe)dlog Eg. (5) 

Since the coronal abundance of low-FIP elements is unknown a priori, we recast Equation ([5]) in 
terms of the photospheric abundances which are believed to be known. Thus, the observed Fe 
intensities are 

aP r+oo 

IZ = 17 y^ ng)ep(Eg)dlog Eg, (6) 

where we use the relationship Apg = /pgA^ and is the coronal DEM derived using photospheric 
abundances. The DEM is related to the coronal DEM by 


ep(Eg) = /FgCc(Eg). 


(7) 


Using the approach of lHannah fc Kontarl (120121 1 we can invert Equation ([^ to find the inferred 
^p. Figure [3] displays a typical example of such a DEM ^p and will be discussed in more detail in 
Section 13.41 


3.3. The Fe-to-Si FlP-bias ratio 

The next step is to compute the FlP-bias ratio /pe/Zsi- Since Si and Fe are both low-FIP 
elements, with similar EIP values, they are expected to have a similar EIP bias. Numerous 








observ ational studies have demonstra t ed this in measuremen ts of many different coronal struc¬ 
tures ([Feldman &: Widina l2003l . 120071 : iBrooks &: WarrenI l201lh . In order to estimate the fpe/fsi 
ratio, we compare the predicted and observed intensities of the 6 Si lines used in this analysis. 


The predicted Si intensities are computed using the photospheric silicon abundance. As before 
the photospheric and coronal abundances are related by the silicon FIP bias fsi as Agj = /siAgj. 
Using the inferred DEM derived from the Fe lines by inverting Equation ([6|), the predicted Si 
intensities can be written as 

r+oo 

Iwed=fFe^J^ G,,{T,,n,)UTe)d log T„ (8) 

where has been recast using Equation ([7j). The observed Si intensities are given by 


r+oo 

= /Si^ G',,( Te, ne)ec( Te)d log Te. 

Therefore, the ratio of the predicted to observed Si intensity reduces to 


7-Si 

pred 

-'obs 


/Pe 

/si 


(9) 


( 10 ) 


If both Ee and Si have the same FIP bias, then the observed and computed Si intensities should 
match. 


The question now arises of which compilation of photospheric abundances to use for our anal¬ 
ysis. Several sets are available. In order to determine which set of the reported photospheric 
abundances is consistent with o ur observations , we t ested the effect on our measu rements us¬ 
ing the recommended values of iGrevesse et al.l ([20071 1 and of lAsolund et al.l ([2009l l . The sili¬ 
con abundance measured in these two studies is identical wh ereas that for i ron d iffers, having 
= 1-122. We find that the abundances of Aspluud et ah ( 2009) give a . FIP - 
bias ratio /pe/Zsi of close to unity. This is as expected, and so we use the lAsplund et al.l ([20091 ) 
photospheric abundances for the rest of the analysis. 


Figure 0] shows the variation of the fpe/fsi ratio for each of our three data sets. The white 
areas correspond to rejected spatial bins, following the criteria we described in Section [3Tl We also 
rejected data for which the radius is i? < 1.02 Rq. At such heights, spicules are present, which can 
lead to a multithermal distribution and therefore confuse the analysis. The corresponding profiles 
as a function of the solar-X coordinate are displayed in Figure (U along with the intensity profiles 
observed using the Fe X 184.54 A line to identify the position of each plumes. For each data set, 
we find a relatively homogeneous /pe/Zsi ratio of between 0.9 and 1.1. There is no clear correlation 
between the measured /pe//si ratio and the Fe X intensity. It is worth noting that these maps are 
sparser than those presented later in this work. This is because we have used only iron lines to 
compute the DEM, and so fewer spatial bins were consistent with the criteria that we applied to 
ensure the reliability of the DEM inversion. 
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3.4. The Fe & Si DEM 

The /pe/Zsi ratio is close to unity, which implies that Fe and Si have a similar FIP bias. This 
motivates our computing a Fe & Si DEM using the Fe and Si lines for each spatial bin. We used 
Equation @ for the Fe lines. For the Si lines, we used a similar equation but with in place of 
Agj. This allows more lines to be used in the inversion process, thereby providing a more robust 
DEM estimation. The DEM, inferred using photospheric abundances, is now given by 

?p(?e) = /(Fe & Si)^c(Fe), (11) 

where /(pe & si) = /pe = /si- Figure [3] displays a typical example of both Fe and Fe & Si DEMs. 
The stars represent the quantity Cp(^t).fpred/.^obs for the Fe, Si, and S lines, computed using the Fe 
&: Si DEM. Tt is the DEM weighted formation temperature, given by 


_ /G,,(re)gp(re)logTedre 
* /G,i(re)^p(Te)dre 


( 12 ) 


The quantity Cp(^t).fpred/.^obs illustrates the match between the inferred DEM and the intensity 
data. 


The Fe and the Fe &: Si DEMs are very similar, though there are some differences. We attribute 
these, in part, to the different number of lines used and thus to the different temperature coverage. 
In general we find good agreement between the predicted and observed intensities, as illustrated in 
Figure [3] by the closeness of most of the stars to the DEM curve. 


The main contribution to the DEM comes from plasma between logTe = 5.9 and 6.1, which 
is the characteristic temperature range for plume and interplume plasma. There is also a small 
contribution at high temperature, between logTe = 6.2 and 6.3, coming from quiet Sun mate¬ 


rial l ying along the LOS. This signal complicates the interpretation of the results (jHahn et al 


20111 ) and can distort the abundance measurements as the quiet Sun is subject to a small FIP 


effect (jFeldman fc Widind 120031 . 120071 ). To avoid this contamination, we developed an additional 
test, based on Pearson’s test, that rejects spatial bins for which this quiet Sun contamination 
is excessive. First, we generated a Fe &: Si DEM model from the DEMs inferred from the first 
data set. For this particular observation, there is little or no contamination from the quiet Sun. A 
DEM model, is computed by averaging the normalized inferred DEMs over all gas at 

Y = —1007'b This average DEM is plotted in Figure [6l Next, for each spatial bin of each data set 
we applied Pearson’s tost in order to evaluate how much the inferred Fe &: Si DEMs differ from 
the model. We then reject Fe & Si DEMs with a significant high temperature contribution. For 
this, we computed Pearson’s cumulative test statistic 


2_^ Mgp(^e')/Cr^-Cm(Ti)]^ 

^ h) UiTi) 


( 13 ) 
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where j indexes a point on the temperature grid and A^bin = 26 is the number of points on this 
grid. If the inferred DEMs are only affected by normally distributed random errors, then Pearson’s 
cumulative test statistic is equivalent to a statistical chi-squared with A^bin degrees of freedom. 
Spatial bins are rejected if is too large. We define the threshold for this to be > 38.885. The 
probability of have such a large x^ for a statistical 26-degrees-of-freedom, chi-squared distribution 
is less than 5%. That implies a 95% chance that there is significant quiet Sun contamination in the 
spatial bin. 

Figure [3 shows the maps for each of the three data sets. The hrst and second maps exhibit 
very homogeneous x^ values, with most of the values falling between 0 and 5, indicating that the 
observations are of relatively pure coronal hole plasma. However, the last observation clearly shows 
that some of the inferred DEMs contains significant quiet sun contamination. These spatial bins 
have been rejected. 


3.5. Computing the Fe & Si to S FlP-bias ratio 


In order to derive the relative FIP bias, we compared the observed S VIII intensities to the 
intensities predicted using the inferred Fe & Si DEM ^p. From this DEM, derived in Section [3.41 
the predicted intensities are given by 

r+co 

-^pred ~/(Fe & Si) J Gji{Te, ne)^p{Te) d log Tq. (14) 

The observed S VIII 198.55 A intensities are given by 


j^P |• + CO 

iL = fs^ Ue)^p(re)dlog Te, 


(15) 


where /s is the sulfur FIP bias. The ratio of the predicted-to-observed S VIII 198.55 A intensities 
is simply 


rS 

pred 



/(Fe fc Si) 

/s 


(16) 


Using this for each spatial bin of our three observations, we have generated maps of /(pe & Si)//s 
which we discuss in Section [H 


4. Results 

Figure [8] shows the /(pg ^ si)//s FlP-bias-ratio maps for the three data sets. The correspond¬ 
ing profiles, computed for the maps given in Figure [8] by averaging the spatial bins over the Y 
direction, are displayed in the Figured Also shown in Figure [9] are the intensity profiles observed 
using the Fe X 184.54 A line. Table [2] lists the identified plume and interplume regions in each 
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observation; the approximate X location of their intensity maximum or minimum, repsectively; 
and their approximate FlP-bias ratios. 

The top panel of Figure [8] shows /(pg ^ Si)//s for the first data set, while the corresponding 
profiles are shown in the top panel of Figure [H Two interplume regions are visible at X ~ —50^^ 
and 60^h Three plumes are easily identified in this particular observation, located respectively at 
X ~ —15'^, If)” and 85'h These coincide with the intensity enhancement seen in the contour lines 
of Figure [8] and the intensity profiles of Figure [9l There is a suggestion of a fourth plume, located 
at X ~ —85" and visible on the FIT field of view of Figure [2l 

The two interplume regions present a FlP-bias ratio of ~ 0.7. Conversely, in the plume regions 
the FlP-bias ratio varies from plume to plume. The plume located at ~ —15" exhibits the largest 
relative FlP-bias ratio, with values close to 0.95. The plume visible at ~ 85", shows an enhancement 
of about 0.8. On the other hand, the brightest plume, detected at « 15", has a FlP-bias ratio of 
about 0.7, similar to what we measure in the interplume regions. We save discussion of the plume 
at ~ —85" until later. 

The middle panels of Figure [8] and [9] display the /(pg & si)//s measured for the second data 
set. Even though the field of view is smaller, we still clearly observe the plume located around 
X R3 —20". The brightest plume, detected at around X « 15" in the top panel, has faded away, 
and an interplume region is now visible at X ~ 10". The plume still exhibits a FlP-bias ratio of 
about 0.9. This is slightly smaller than in the first data set, but still larger than the FlP-bias ratio 
of the interplume region, which is about 0.7, similar to the results from the first data set. 

The bottom panels of Figure [8] and [9] present the last data set, collected 24 hours after the 
first. Two plumes are still clearly visible, now located at X ~ —30" and 100". We also find a 
suggestion of a third plume at X ~ —100", which is rather fainter than the other two. This may be 
the plume from the first data set at X ~ —85". We attribute the shifts in the positions of all the 
observed structures to a combination of the solar rotation and movement of the held lines caused 
by other possible processes. Interplume regions are detected around the X —70" and 50", and 
continue to exhibit FlP-bias ratios of about 0.7. The plume around X ~ —30", which displayed 
the largest enhancement in the hrst two data sets, now exhibits a reduced FlP-bias ratio of about 
0.8. This is particularly clear when looking at the prohles in Figure [9l Also notable is that all three 
of the remaining plumes have a similar FlP-bias ratio of about 0.8. 

To summarize this section, we hnd that interplume gas exhibits a constant FlP-bias ratio 
over the duration of the observations. Conversely, over the ~ 24 hour observation period the time 
evolution of the FlP-bias ratio differs dramatically from plume to plume. The initially brightest 
plume at X ~ 15", which appears to originate from a bright point, rapidly fades away, becoming an 
interplume region, and does not show any evolution in its FlP-bias ratio, with an initial FlP-bias 
ratio similar to that of the interplume regions. This is to be contrasted with the behavior of the 
other three plumes. The plume initially at X ~ —15" shows the largest initial FlP-bias ratio, but 
the value decreases over time, approaching, though not reaching, that of the interplume gas. These 
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changes occur despite the apparently nearly constant brightness and density of this structure. The 
plume initially at X ~ —80" shows a small increase in its FlP-bias ratio over time, increasing from 
a value initially close to that of interplume gas. Lastly, the plume initially at X ~ 85" exhibits a 
nearly constant FlP-bias ratio, which is slightly higher than that of the interplume gas. 


5. Systematic Uncertainties 

There are several source of uncertainties involved in our analysis. In this section we attempt 
to identify these sources and discuss their possible effects on our results. 


5.1. Instrumental 


Instrumental factors, such as the ca libration of the EIS instrument could produce a systematic 
bias in the results. ICulhane et al.l (j2007l l estimated the absolute calibration of the EIS spectrometer 
to be around 25% for each of the CCD cameras on the instrument. However the relative calibration 
error within a given camera or between the two cameras is expected to be much smaller. Thus 
we expect this error essentially cancels out in our analysis. Additional instrumental issues could 
arise due to degradation effects. However, we used data from 2007, at the be ginning of the EI S 
observations period, when there had not yet been any significant degradation (iDel Zannall201.‘ll ). 
Thus no significant instrumental uncertainties are expected in our results. 


5.2. Atomic data 


Errors associated with the atomic data required to interpret the observed intensi ties are another 


poten tial source of systematic uncertainty, but such errors are difficult to estimate ([Eeldman et al 


20091 ) . One source of error could be in the line emissivities due to uncertaintie s in the exci t atioii rate 
coefficients, transition rates, and incompleteness of the cascade contributions. iDel Zannal (j2018l ) has 
carried out a detailed study of solar observations compared to model predictions and identified a 
number of lines where the model line emissivities are likely to be in error. We have used their 
results as a guide and avoided all those lines. Hence, we expect the line emissivity uncertainties to 
be minimal for the lines used in the present study. 

Another source of atomic physics uncertainty lies with the ionization balance calculations. 
T here has been e xtens ive experimental and theoretical work to improve these data as summarized 


m 


Brvans et al.l (120091 ). Our results her e used the ionizatio n balance calculations of CHIANTI 


version 7.1 which has updated the data of iBrvans et al.l (j2009l ) by including state-of-the-art recom¬ 
bination rate coefficients for isoelectronic sequences published after their work. Still, in order to 
estimate the effects of the new ionization balance calculations on our measurements, we repeated 
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the entire analysis using the data from iBrvans et alJ (|2009l l. We did not find any substantial change 
from the results presented here. The difference between the two FlP-bias-ratio maps never exceeds 
5%, and the relative variations of the FlP-bias ratio between plumes and interplumes is unchanged. 


5.3. Intensity ratios: the isothermal approximation 


The intrinsic differences between the plume and interplume DEMs might introduce a systematic 
bias in the inferred FlP-bias ratio. Figure [TU] displays typical inferred Fe &: Si DEMs for plume 
and interplume gas. Both functions peak around logTg = 5.95, but the interplume peak is weaker 
than the plume, due to its lower density. The plume DEM also exhibits a higher contribution from 
plasma between logTe = 5.8 — 5.9, coming from the cooler part of the plume plasma. 

We hnd that the DEMs measured in both plumes and interplumes are narrow in temperature. 
Within the Gaussian approximation, the full width at half maximum (FWHM) is AlogTg ~ 0.15 
for the interplume regions, and A logTe ~ 0.2 in plume regions. T hese are small value s , con ¬ 
side ring the temperature r esolution of the DEM inversion. Indeed, iLandi fc KlimchukI ()201fll ) 
and iGuennou et al.l (l2012al ) have demonstrated that plasmas with a Gaussian-DEM FWHM of 
A logTe < 0.10 are consistent with the gas being isothermal. This is due to the presence of sys¬ 
tematic and random errors, which limit the intrinsic temperature resolution of the DEM inversion. 


The narrow widths of the DEMs for both the plumes and interplumes enables us to use the 
isothermal approximation in order to estimate the error on the absolute ElP-bias-ratio results. This 
also allows us to test the robustness of the relative variation in the inferred FlP-bias ratios. For 
this we use the Si VII 275.39 A and S VIII 198.55 A lines, which are particularly well suited for 
FlP-bias-ratio measurements of coronal hole plasma, since most of their line emission is from the 
temperature ranges typical of coronal holes (jFeldman et al.M2009l ). 


In the isothermal approximation, the observed intensity of the Si VII 275.39 A line for a plasma 
at a given temperature Te and density Ue is 


-^obs =/si 




a (^si^ 


dvr 


p 


Te, Ue) X EM, 
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where G®) is the Si contribution function and EM is the emission measure, which is defined using 


Equation ([2]) as 


I 


EM = / nidi. 


In the same way, the S VIII 198.55 A intensity can be written 

4s=/s^G'|(Te,ne)xEM, 


(18) 
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where Gf,- is the S contribution function. 

J ^ 
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The intensity ratio of the Si VII 275.39 A to the S VIII 198.55 A lines gives 

4, /s.4g Gf;(r„n,) 

4. /siAg,0|i(r.,n.)’ 

and thus the FlP-bias ratio /si//s can be calculated using 

A, 

fs .4P Gj.(T,.ne)4.' 

In order to estimate the FlP-bias ratio using Equation (|2ip . we need to calculate the contribution 
function for the Si VII 275.39 A and S VIII 198.55 A lines. For this, the density has been estimated 
using the Fe IX line intensity ratio as described in Section 12.31 Additionally, for each spatial bin, 
we estimated Tg as the temperature corresponding to the peak of the DEM. 

Figures [11] and [12] respectively show the resulting FlP-bias-ratio maps and the prohles of the 
FlP-bias ratio fsi/fs for our three data sets. Table [3] lists the identified plume and interplume 
regions in each observation; the approximate X location of their intensity maximum or minimum, 
respectively; and their approximate FlP-bias ratios. 

The top panels of Figures [TT] and [12] show that the two interplumes regions, located at A ~ 
—50" and 60", have FlP-bias ratios of 1.05 and 0.97, respectively. The two plumes located at 
X ~ —15" and 85" exhibit higher FlP-bias ratios of ~ 1.32 and 1.12, respectively. Conversely, the 
plumes at A ~ —80" and 15" have FlP-bias ratios slightly smaller than those of the interplumes 
regions. In general, these relative results are similar to those of the DEM analysis, though the 
absolute values of the ElP-bias ratios are different. 

In the middle panels, corresponding to the second data set, the plume is now located at 
A ~ —20" and shows a ElP-bias ratio of ~ 1.23. The newly formed interplume region at A ~ 10" 
exhibits a FlP-bias ratio ~ 1.12, which is actually higher than the now-faded-plume results of the 
first data set. Thought the differences between the two regions are less marked than in the DEM 
analysis, the general trend is similar. 

The bottom panels of Figures [TT] and W2\ correspond to the third data set. The FlP-bias ratio 
of the plume located at A ~ —100" has increased slightly, similar to what was seen in the DEM 
analysis; however, the ratio is not significantly different from that of the interplume results. For 
the plume at A ~ —30", we find a decreasing trend in the FlP-bias ratio over the course of the 
three observations, similar to what was seen in the DEM analysis. Lastly, the FlP-bias ratios of 
the plume located at A ~ 100" is about 1.25, higher than that of the interplume gas in the first 
observation. However, this plume exhibits a different behavior than was inferred from using the 
DEM analysis. In the isothermal case the FlP-bias ratio increases, whereas in the DEM analysis it 
remained nearly constant. 

The isothermal analysis, similar to that of the DEM, finds that interplume gas exhibits a 
roughly constant FlP-bias ratio over the duration of the observations. Additionally, over the ~ 24 
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hour observation period the time evolution of the FlP-bias ratio varies from plume to plume. De¬ 
spite the minor differences between the isothermal and DEM analyses of the plumes, this strengthens 
the conclusions drawn from the DEM results. However, we still need to address the uncertainties 
in the inferred FlP-bias ratios. 

We used a Monte-Carlo procedure to estimate the uncertainty associated with the isother¬ 
mal measurements of the FlP-bias ratio fsi/fs- For this we took into account the errors on the 
intensity, density, and temperature measurements. In order to estimate the effect from each of 
these uncertainties, we randomly varied one of the parameters, using a normal distribution and 
the associated la error, while keeping the other two parameters fixed. The modified parameter 
is used to estimate the FlP-bias ratio, and the operation is repeated 500 times. The impact of 
the uncertainty in each parameter on the FlP-bias ratio is deduced by computing the standard 
deviation of these 500 different results. The individual uncertainties associated with the intensity, 
density, and temperature measurement errors are reported in Table [3j 


5.3.1. Uncertainty due to the Intensity and Density Measurements 

The statistical uncertainties in the intensities have been computed as described in Section [2.21 
There may also be some systematic uncertainties related to the absolute calibration of the EIS spec¬ 
trometer (see Section [5T]) . Since the FlP-bias-ratio calculation takes the ratio of the Si VII 275.39 A 
and S VIII 198.55 A line intensities, this systematic error essentially cancels out and only the small 
relative calibration error remains. The impact of this relative error on the FlP-bias-ratio estimation 
is expected to be insignificant. 

We find that the statistical errors involved in the intensity measurements give rise to a Icr 
uncertainty in the FlP-bias ratio of between 0.08 and 0.12, varying with the location in the obser¬ 
vations. In general, this is smaller than the observed FlP-bias-ratio differences. Hence we conclude 
that uncertainties in the intensity measurements cannot explain the observed structures. 

The uncertainty in the density has been estimated using the computed Icr errors on the 
Fe IX 188.50 A and 189.94 A lines intensities. As described in Section 12.31 the density diag¬ 
nostic assumes that the plasma is isothermal at the temperature corresponding to the maximum 
of the Fe IX 188.50 A and 189.94 A contribution functions, which peak at about logTe ~ 5.9. 
This temperature is essentially indistinguishable from the peak temperature of the DEMs for both 
plumes and interplumes. Hence, we expect any small systematic errors due to the temperature 
used to be insignificant. Additional systematic uncertainties could arise due to errors with the ex¬ 
citation rate coefficients, transition rates, and incompleteness of the modeled cascade contributions 
to the Ee IX 188.50 A and 189.94 A line intensities. Fortunately, the contribution functions for 
the Si VII 275.39 A and S VIII 198.55 A are only slightly dependent on the density. As a result, 
systematic uncertainties related to the density measurements are expected to only have a small 
impact on the FlP-bias ratio calculations. We therefore only consider the statistical error related 
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to the intensity measurements. 

We find that errors in the density measurements give rise to a la uncertainty in the FlP-bias 
ratio of less than 0.02 for both plume and interplume gas. This is smaller than most of the observed 
FlP-bias-ratio differences. Hence we conclude that uncertainties in the density measurements can¬ 
not explain the observed structures. 

The quadrature sum of the uncertainties due to the intensity and density measurements gives 
a Icr value of less than 0.12. Again, this is still smaller than most of the observed FlP-bias-ratio 
differences and strengthens our conclusion that the observed structure cannot be attributed to 
errors in the intensity or density measurements. 


5.3.2. Uncertainties due to the Temperature Measurements 


The temperature uncertainty is taken from the work of iLandi &: KlimchukI (|2010l l and iGuennou et al. 


(|2012al l. They showed that, due to DEM-inversion complications, the uncertainty of the inferred 
DEM-peak temperature is estimated to have a EWHM of A log Eg = 0.1. This takes into account 
both the random and systematic unc ertainties due to unknown ato mic physics and i nstrum ental 
calibration. However, as described in iLandi &: KlimchukI (|2010l i and iGuennou et al.l (|2012al l. the 
contribution of the statistical uncertainties to this estimated temperature error is much smaller 
than that due to systematic causes. 


The FlP-bias uncertainty related to the temperature error is higher than that due to the 
intensity or density. The Icr value ranges between 0.46 and 0.57. This is as expected, since the 
contribution functions are strongly dependent on the temperature. 


At first glance, it appears that this uncertainty is large enough to wash out the structure seen 
in the three data sets. However, we find it unlikely that the temperature in plumes and interplumes 
would be mis-estimated in such a way as to artificially generate the observed structure. If such 
a systematic effect were at work, then all of the plumes would be likely to show a ElP-bias ratio 
different from that of the interplumes. What we see instead is that only some of the plumes differ 
from the interplumes in their ElP-bias ratio. Additionally, the DEM analysis effectively takes into 
account the temperature uncertainty in the isothermal analysis; but it shows a structure similar to 
the isothermal results. Hence we interpret the temperature uncertainty in the isothermal results to 
represent a systematic uncertainty in the absolute FlP-bias results. We also take it as an estimate 
for the systematic uncertainty in the DEM ElP-bias results. 


5.3.3. Robustness of the Results 

Based on the above error analysis, we have a high degree of confidence in the robustness of our 
relative FlP-bias-ratio results. The uncertainties due to the intensity and density measurements are 
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too small to explain the observed relative differences; and the uncertainty due to the temperature 
is most likely to systematically affect the results. Taking this last uncertainty as an estimate for 
the systematic reliability of the absolute scale, we find that this ~ ±0.5 accuracy of the absolute 
scale can readily explain the difference seen between the DEM FlP-bias ratio /(pg ^ Si)/fs, which 
ranges between ~ 0.7 and 0.9, and the isothermal FlP-bias /si//s, which ranges between ~ 1 and 
1.3. 


6. Discussion 


Over the 24 hour period studied, in the interplume regions the DEM-inferred FlP-bias ratio 
/(Fe & Si)//s exhibit a nearly constant behavior. Conversely, the plume regions display a range of 
behaviors. One starts with a nearly interplume FlP-bias ratio, which then increases slightly with 
time. The plume with the highest initial FlP-bias ratio shows a decreasing ratio with time, though 
it never reaches interplume values. A third plume begins with an interplume FlP-bias ratio and 
then fades, becoming an interplume region with nearly the same initial FlP-bias ratio. The fourth 
plume shows a nearly constant FlP-bias ratio, but one which is great than the interplume value. We 
have a high degree of confidence in the robustness of these relative results. Our DEM findings are 
largely supported by our corresponding isothermal analysis. Thus, we conclude that the observed 
variations are most likely due to real compositional differences. 


Building on our findings , the evolution of plumes with time could explain the discrepancie s 


seen among previous studies (jWilhelm fc Bodmer 


1998; 


Young et al.lll999l : iDel Zanna et al. 


200311 . 


These past measurements may have observed plumes at different stages in their lifetime. Although 
only one plume in our study showed a clear elemental abundance variation over the ~ 24 hr period, 
our findings demonstrate that abundance evolution can occur in plumes. Clearly, though, the others 
remained essentially stable in time. More systematic studies are needed in order to determine what 
are the common compositional features of plumes. 


One possible explanation for the observed time evolution of the FlP-bias ratio is that it is 
related to the time at which the plume formed over the underlying bright point. Beam plumes 
are believed to originate from magnetic reconnection of em erging magnetic flux with the am bient 
unipol ar magnetic fie l d, the reby generating a bright point (jWandll994l : iGrappin et al.l 120111 ). Re¬ 
cently, iRaouafi et al.l (120081 ) showed that jets that arise from magnetic reconnection are precursors 
of plumes. Previous studies have shown that the FIP bias is related to the confinement time of the 
plasma. Thus, the initial abundance enhancement in the plume may represent the FIP bias in the 
bright point at the time the jet occurred. Once the material confined in the bright point is ejected, 
photospheric material is then directly ejected through the open magnetic field lines, resulting in a 
subsequent change of the elemental abundances towards photospheric, as observed in our results. If 
this hypothesis is correct, FlP-bias ratio measurements, such as those presented here, can provide 
constraints on the confinement time of the plasma that is released when the plumes are formed. 
Such data may also provide important boundary conditions for models of the plumes formation 
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process and the heating mechanism occurring at their base. 

A possible future application of the approach laid out in this work is to be able to identify 
plume and interplume material in the solar wind. Comparing the present results to the available 
solar wind in situ measurements is complicated, because our FlP-bias-ratio measurements have 
significant systematic uncertainties. However, archival spectral observations from SOHO/CDS are 
available. These data include emission lines from many more elements and ionization stages than 
do the EIS data. The CDS data would thereby enable the construction of better constrained 
DEMs, resulting in reduced uncertainties in the inferred ElP-bias-ratio measurements. We hope in 
a future work to use CDS data for such studies, since more plume studies are needed in order to 
better characterize any time evolution of elemental abundances and develop diagnostics so as to be 
able to identify the signature of plumes within the fast solar wind. 


At the moment, the Sun is in its active period, with minimally developed coronal holes and 
excessive quiet Sun contamination, thereby making present-day plume and interplume observations 
more difficult. However, future observations from spectrometers and in situ instruments on board 
Solar Orbiter offer great opportunities for connecting polar plume and interplume observations to 
in situ solar wind measurements. Such observations are likely to help to answer the question of 
whether or not plumes contribute to the fast solar wind. 


Lastly, our DEM analysis of the plume and interplume regions reveal a ElP-bias ratio /(pe & Si) / /s 
that is systematically less than 1.0. This would imply either a depletion of the low-EIP elements 
or an enhancement of the sulfnr. A depletion w ould suggest an inve r se EIP effect, which t o date 
has only been observed in the most active stars ( Feldman et al. 2001 : Brinkman et aklbooih . Pre¬ 
vious FIP bia s measurements of coronal hole plasma have consistently found a low-FIP element 
enhancement ([Feldman et al.ll2005l : iFeldman k, Widindl2007l : ICrevesse et al.ll2007ll . As discussed in 
Section [5l we attribute the ratio here of less than 1 to a systematic error due to the uncertainty in 
the DEM shape. Euture measurements could possibly reduce this error by using more lines covering 
a wider temperature range, thereby enabling a more accurate determination of the DEM. 


7. Summary 

In this work, we have described measurements of the ElP-bias ratio /(pe & si)//s within plume 
and interplume plasma. We used spectroscopic observations over an ~ 24 hour period in order 
to characterize the evolution of these structures with time. We have found that the ElP-bias 
ratio is essentially constant with time in interplumes but can vary with time in plumes. Euture 
studies combined with corresponding in situ measurements of the solar wind may help to answer 
the question of whether or not plumes or interplumes contribute to the fast solar wind and may 
also provide constraints on possible formation and heating mechanisms of plumes. 
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Table 1. List of spectral lines used. 


Ions 

Wavelength (A) 

iog(r[iL]) 

Si VII 

272.65 

5.70 

Si VII 

275.36 

5.70 

Si VII 

275.67 

5.70 

Si VII {bl Mg VII) 

278.45 

5.70 

Six 

261.06 

6.10 

Six 

277.26 

6.10 

Fe VIII 

185.21 

5.70 

Fe VIII 

186.60 

5.70 

Fe VIII 

194.66 

5.70 

Fe IX 

188.50 

5.75 

Fe IX 

189.94 

5.75 

Fe IX 

197.86 

5.75 

Fe X 

184.54 

6.00 

Fe X 

190.04 

6.00 

Fe XI 

182.17 

6.00 

Fe XI 

188.22 

6.00 

Fe XI 

189.12 

6.00 

Fe XII 

192.40 

6.10 

Fe XII (sbl) 

195.12 

6.10 

Fe XII 

196.54 

6.10 

Fe XIII 

202.04 

6.15 

Fe XV 

284.16 

6.30 

S VIII {bl Fe XI) 

198.55 

5.80 


Note. — Emission lines are sorted by elements 
and ordered as a function of the peak temperature 
of their contribution functions. The blended and 
self-blended lines are indicated by the abbreviations 
bl and sbl, respectively. 
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Table 2. Identified plumes and interplumes; the location of their intensity maximum or 
minimum, respectively; and their FlP-bias ratio, /(pe & Si)//S) for the DEM analysis. 


Data set 



1 


2 


3 


Region Location 

/(Fe & Si)//S 

Location 

/(Fe & Si)//S 

Location 

k 

Fe & Si)//s 

Plume 

1 

00 

o 

^ 0.73 



o 

o 

r—1 

1 


0.78 

Interplume 

« -50" 

^ 0.71 



Rs -70" 


0.71 

Plume 

« -15" 

ss 0.94 

o 

1 

« 0.88 

o 

CO 

1 


0.78 

Interplume 



« 10" 

« 0.72 

5" 


^ 0.71 

Plume 

15" 

RS 0.71 






Interplume 

^ 60" 

0.69 



50" 


0.70 

Plume 

85" 

ss 0.80 



« 100" 


0.79 






Table 3. Identified plumes and interplumes; the location of their intensity maximum or minimum, respectively; and their 
FlP-bias ratio, fsi/fs^ follow by the la errors for the intensity, density and temperature, respectively, for the isothermal 

analysis. 






Data set 





1 


2 


3 

Region 

Location 

/si//s 

Location 

fsi/fs 

Location 

fsi/fs 

Plume 

1 

00 

o 

0.95 ±0.08,0.01,0.46 



« -100" 

1.08 ±0.08,0.01,0.49 

Interplume 

« -50" 

1.05 ±0.10,0.01,0.54 



^ -70" 

0.96 ±0.11,0.02,0.52 

Plume 

« -15" 

1.32 ±0.10,0.01,0.60 

o 

1 

1.23 ±0.08,0.03,0.53 

O 

CO 

1 

1.09 ±0.12,0.02,0.61 

Interplume 

Plume 

^ 15" 

1.02 ±0.09,0.01,0.51 

« 10" 

1.12 ±0.08,0.02,0.46 

^5" 

0.99 ±0.09,0.01,0.52 

Interplume 

^ 60" 

0.97 ±0.11,0.01,0.50 



«50" 

1.03 ±0.10,0.02,0.46 

Plume 

^85" 

1.12 ±0.10,0.02,0.57 



« 100" 

1.25 ±0.11,0.02,0.53 
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Fig. 1.— The box outlines the raster range covered by the EIS V slit observations on 2007 March 
13 at 13:00 UT. The solar image is from the S'0770/FIT 171 A observation closest in time to the 
EIS raster and consists mainly of emission from Fe X and Fe XL 
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Fig. 2.— Left column: Maps of the measured intensity for each of the three data sets as seen by 
50if0/EIT in the 171 A channel. The red box on the middle pannel outlines the position of the 
EIS raster for this data set. Right column: Maps of the measured density on a logarithmic scale 
for each of the three EIS data sets. The black contours correspond to the intensity measured using 
the Fe X 184.54 A line. The solar limb is highlighted by the white line. 
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Fig. 3.— Example of a typical .^p DEMs computed for one spatial bin at X ~ 76" and Y ~ 
—997'', corresponding to an interplume area. Both the iron DEM (orange line) and the Fe & 
Si DEM (blue line) are displayed. For this particular spatial bin, 13 Fe and 3 Si lines match 
the criteria described in Section 13.11 ensuring that the DEM is well constrained. The horizontal 
error bars represent the temperature resolution of the inve rsion and the vert i cal er ror bars are 
the uncertainties arising from the intensity measurements (iHannah fc Kontarl 120121 ) . Each star 
represents the quantity ^p{Tt)Ipred/Iohs computed for the Fe & Si DEM (see Section [TH for details). 
There is good agreement between the Fe and the Fe & Si DEMs, indicating that both Fe and Si have 
a similar FIP bias. Also notable is the very good agreement between Fe and Si in the temperature 
range where S VIII forms (green star). 
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Fig. 4.— Maps of the FlP-bias ratio /pe/Zsi for each data set. The black contours are as in Figure[2j 
The red line displays the solar limb. The FlP-bias ratio is concentrated between 0.9 and 1.1, with 
mean values and Icr-standard deviations of 0.97 ± 0.12, 0.96 ± 0.09, and 0.96 ± 0.10, going from the 
top map to the bottom. 
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Fig. 5.— FlP-bias ratio profiles /pe/Zsi for each observation (red lines and left axes). Each profile 
has been computed by averaging the spatial bins of the maps in Figure 0] along the solar-E axis. In 
order to simultaneously present the plume positions, the observed intensities in the Fe X 184.54 A 
line are also displayed on the same plot (blue lines and right axes). These intensity profiles have 
been obtained by averaging 7 spatial bins from V ~ —1002" to T ~ —1030" in order to enhance 
the signal and increase the visibility of the features. 





~ 30 - 



log [n (K)] 


Fig. 6.— The normalized Fe & Si DEM model (red line) computed for the first data set, where 
there is little contamination from high temperature emission. The blue line shows an example of 
an inferred, normalized DEM close to the model (x^ = 2.30), located at X ~ —36 and Y ~ —1032. 
The green line shows an inferred normalized DEM with a significant contribution from quiet Sun 
plasma, leading to a x^ = 39.54, for the spatial bin at X ~ 96 and Y ~ —1001. According to the 
X^ test described in Section [3~il we rejected this spatial bin from our analysis. 
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Fig. 7.— Maps of Pearson’s cumulative test statistic for each data set, in the same order as in 
Figure [2l The black regions correspond to spatial bins where x^ > 38.885 and which have been 
rejected due to excessive quiet Sun material. The last observation is the most contaminated by 
high temperature emission, as can be seen on the right hand side of the map. The black contours 
are as in Figure [2] and the red curve as in Figure [H 
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Fig. 8.— Same as Figure 01 but for the FlP-bias ratio /(pe & Si)//s- 
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Fig. 9.— Same as Figure O but for the FlP-bias ratio /(pg ^ si)//s- 
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Fig. 10.— Typical plume (red line) and interplume (blue line) Fe &: Si DEMs. The plume 
DEM is computed for the spatial bin located at X = —8", while the interplume is at X = 57". 
The solar-T coordinate in the both cases is T = —997'h The horizontal error bars represent the 
temperature resolution of the inversion and the vertical error bars are the uncertainties arising from 
the intensity measurements. 
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Fig. 11. — Same as Figures S] and [8l but for the FlP-bias ratio /si//s using an isothermal analysis. 
These maps have been obtained using the line pairs diagnostic Si VII 275.39 A and S VIII 198.55 A 
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Fig. 12.— Same as Figures [5] and [U but for the FlP-bias ratio /si//s using an isothermal analysis. 
These curves have been obtained using the line pairs diagnostic Si VII 275.39 A and S VIII 198.55 A 





